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Abstract 

Background: Pliospliorus (P) deficiency is one of tlie major nutrient stresses limiting plant growth. The uptake of 
P by plants is well considered to be mediated by a number of high-affinity phosphate (Pi) transporters belonging 
to the Phtl family. Although the Phtl genes have been extensively identified in several plant species, there is a lack 
of systematic analysis of the Phtl gene family in any solanaceous species thus far. 

Results: Here, we report the genome-wide analysis, phylogenetic evolution and expression patterns of the Phtl genes 
in tomato {Solanum lycopersicum). A total of eight putative Phtl genes {LePTl to 8), distributed on three chromosomes 
(3, 6 and 9), were identified through extensive searches of the released tomato genome sequence database. Chromosomal 
organization and phylogenetic tree analysis suggested that the six Phtl paralogues, LePTl/3, LePT2/6 and LePT4/5, which 
were assigned into three pairs with very close physical distance, were produced from recent tandem duplication events 
that occurred after Solanaceae splitting with other dicot families. Expression analysis of these Phtl members revealed that 
except LePT8, of which the transcript was undetectable in all tissues, the other seven paralogues showed differential but 
partial-overlapping expression patterns. LePTl and LePT/were ubiquitously expressed in all tissues examined, and their 
transcripts were induced abundantly in response to Pi starvation; LePT2 and LePT6, the two paralogues harboring identical 
coding sequence, were predominantly expressed in Pi-deficient roots; LePT3, LePT4 and LePT5 were strongly activated in the 
roots colonized by arbuscular mycorrhizal fungi under low-P, but not high-P condition. Histochemical analysis revealed that 
a 1250-bp LePT3 promoter fragment and a 471 -bp LePT5 promoter fragment containing the two elements, MYCS and 
PI BS, were sufficient to direct the GUS reporter expression in mycorrhizal roots and were limited to distinct cells harboring 
AM fungal structures. Additionally, the four paralogues, LePTl, LePT2, LePT6 and LePT7, were very significantly down-regulated 
in the mycorrhizal roots under low Pi supply condition. 

Conclusions: The results obtained from this study provide new insights into the evolutionary expansion, functional 
divergence and genetic redundancy of the Phtl genes in response to Pi deficiency and mycorrhizal symbiosis in tomato. 
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Background 

Phosphorus (P) is one of the three most essential macronu- 
trients required by plants. It is well recognized as serving a 
wide range of structural and biological roles, such as 
energy metabolism, signal transduction, biosynthesis of 
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macromolecules, modulation of respiration, photosynthesis 
and other metabolic processes [1]. The primary source for P 
uptake by plants is orthophosphate (Pi) in soil. Due to the 
slow diffusion rate and chemical fixation, P is widely consid- 
ered to be one of the most difficult nutrients for plants to 
forage, and often a major limiting factor to crop yields [2,3]. 

The Pi concentration in soil solution is commonly no 
more than 10 (iM, whereas plant cells need to maintain 
their cytoplasmic Pi concentrations at a millimolar range 
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[4,5], which determines the requirement of metabolic 
energy and specific transport systems for plants to acquire 
Pi from soils [6,7]. In the past decades, aided by systematic 
studies of molecular biology and functional genomics of 
model plants, a great deal of knowledge about the mecha- 
nisms of Pi transport by plants has been accumulated, re- 
vealing that the uptake and subsequently redistribution of 
Pi within plants are mediated by a number of phosphate 
(Pi) transporters with different affinities that located in the 
plasma or organelle membranes [8]. 

The first gene encoding plant Pi transporter (AtPTl) 
was isolated from Arabidopsis [9] and showed high se- 
quence identity to the genes encoding high-affinity Pi 
transporters in Saccharomyces cerevisiae (Pho84) [10] and 
in Glomus versiforme (GvPT) [11]. The later studies fur- 
ther led to the isolation of other eight homologues exhibit- 
ing substantial identities to the AtPTl in the Arabidopsis 
genome [12], suggesting the expansion of Pi transporter 
genes in higher plants during evolution. By now, with the 
completion of whole genome analysis of model plants, 
such as Arabidopsis and rice, dozens of homologous genes 
encoding different affinities and groups of Pi transporters 
have been identified in various plant species by compara- 
tive genomic approaches [13]. Studies on the protein se- 
quences and phylogenetic relatedness revealed that most 
of the Pi transporters identified so far are typical of H^/Pi 
symporters, and could be grouped into the high-affinity 
Phtl family included in the super facilitator superfamily 
(MPS) [14-16]. 

Earlier studies on the regulation and tissue/cellular 
distribution indicated that members of the Phtl family 
in many species are divergent in function and differen- 
tially expressed during plant development or in response 
to different P status [17,18]. The relatively high levels of 
transcripts or proteins of some Phtl genes in roots, es- 
pecially in root epidermis and root hairs, in response to 
Pi deficiency well support a role of these genes in Pi cap- 
ture and uptake [19,20]. For example, in Arabidopsis^ 
eight of the nine Phtl genes were expressed in roots and 
two members, AtPTl and AtPT4, had the highest ex- 
pression levels in response to Pi deficiency. Knock out of 
either of the two genes showed significant defects in P up- 
take under a low Pi supply condition [21,22]. In some 
cases, the transcripts of some Phtl members are more 
widely distributed throughout plant tissues and showed 
less responses to Pi deficiency, providing strong evidence 
to support that some of the Phtl members may be impli- 
cated in the internal mobilization of Pi, such as loading or 
unloading from the xylem or phloem and deposition into 
seeds or other storage organs [19,23-25]. In addition to 
the Pi-responsive Phtl genes, an increasing number of 
arbuscular mycorrhiza-induced Pi transporters belonging 
to the Phtl family have been identified from several plant 
families, and their functions have been repeatedly docu- 



mented to be associated with Pi uptake at the intraradical 
symbiotic interface [26-33]. 

Tomato, a member of the Solanaceae, is not only a 
world-wide major vegetable crop plant, but also a model 
plant for biological and genetic researches based on its rela- 
tively low-copy DNA sequence and the nearly complete 
genome sequencing [34]. Although previous studies have 
characterized the potential roles of a few individual Phtl 
members in tomato [30,35], there is a lack of genome-wide 
analysis of the Phtl gene family in tomato and also in any 
other solanaceous species thus far. Moreover, compared to the 
other model species, such as Arabidopsis from Brassicaceae 
and rice from Gramineae, the evolutionary mechanisms, tran- 
scriptional regulation and possible frinctions of solanaceous 
Phtl genes in Pi acquisition and mobilization still needs to be 
well explored [36-39]. 

In the current work, we reported the genome-wide 
identification and comparative characterization of Phtl 
family genes in tomato and potato, and further investi- 
gated the expression patterns of tomato Phtl genes in 
response to AM fungi inoculation under low- and high- 
P supply condition. The analysis in this study mainly fo- 
cused on the chromosomal organization, phylogenetic 
evolution, tissue-specific expression and regulation of 
each member of the tomato Phtl family. The results ob- 
tained from this study would not only strengthen our 
understanding on the molecular mechanisms underlying 
the evolutionary expansion, conservation and functional 
divergence of the Phtl genes in tomato, but also provide 
valuable clues for the further comparative genomic stud- 
ies across the whole Solanaceae family. 

Results 

Identification of Phtl family genes in tomato 

Previously, five Phtl genes (three with full-length and 
two with partial mRNA sequence) encoding for putative 
high-affinity Phosphate (Pi) transporters (PT) in tomato 
have been reported [30,40]. In order to determine 
whether there are any further members, as yet unidenti- 
fied, comprising the tomato Phtl family, the mRNA and 
amino acid sequences of Arabidopsis and rice Phtl genes 
were employed for BLASTN and TBLASTN searches 
against the recently released tomato genomic sequence 
database (http://solgenomics.net/), which resulted in the 
identification of a total of eight non-allelic sequences as 
the putative tomato Phtl genes (Additional file 1). BLAST 
searches of these sequences against the NCBI database 
demonstrated that five of the eight sequences were identi- 
cal to the accessioned tomato Phtl genes, LePTl to 5. The 
rest three putative genes (named as LePT6 to 8), which 
were newly identified in this study, showed high levels of 
sequence identity to the known Phtl genes from tomato 
and other plant species (Table 1). Moreover, LePT6, which 
represents a distinct locus, harbors its coding sequence 
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Table 1 Identity matrix for the eight putative tomato Pht1 genes 



Amino acid identity (%) 
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LePT5 


60.9 
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60.1 
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LePT6 


71.0 
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59.3 


57.6 




73.9 


73.1 




LePT7 


74.6 


67.1 


70.5 


59.0 


59.9 


67.1 




92.5 




LePTS 


75.2 


67.0 


69.5 


58.1 


59.0 


67.0 


91.4 





identical to that of the l<nown LePT2, but with much dif- 
ference in un-translated regions between the two homo- 
logues (Additional file 2). 

Comparative analysis of the full-length deduced poly- 
peptides revealed that the eight Phtl proteins contain 
528-538 amino acids with 12 predicted transmembrane- 
spinning domains, similar to the molecular feature of 
Phtl transporters from other plant species. Additionally, 
all the tomato Phtl amino acid sequences share the 
consensus sites for phosphorylation by protein kinase C 
and casein kinase II and conserved residue for N- 
glycosylation (Figure 1). Using the DNAMAN multiple 
sequence alignment program, the conserved domain, 
GGDYPLSATIxSE, which have been suggested to be a 
typical signature of Phtl proteins, was also identified in 
all of these proteins (Figure 1). These findings led to the 
suggestion that all the identified genes could be consid- 
ered as tomato Phtl genes. 

A further blast searches against the tomato EST data- 
base at NCBI, SGN and TIGR database revealed that ex- 
cept LePTSy the other two newly identified Phtl genes, 
LePT6 and LePT7, could matched perfectly to at least one 
significant EST sequences, indicating that the two genes, 
like their previously reported five paralogues, are tran- 
scriptionally active in a certain tissues. It should also be 
emphasized that except the eight Phtl genes mentioned 
above, another sequence (named as LePTx in this study) 
identified in the tomato scaffold database searches also 
showed very high identity to the tomato LePT7 and LePTS 
genes, but may be inactive due to the inclusion of some 
nonsense mutations and indels (insertions and deletions) 
within its putative coding region (Additional file 3), as well 
as to the absence of any EST sequence exactly matching. 

Identification of tomato Phtl homologues in potato and 
comparative analysis of these genes between the two 
solanaceous genomes 

As a near complete set of potato genome sequences 
were also recently available at the SGN database [41], 
for further investigating the evolutionary conservation 



and divergence of Phtl gene family between the two so- 
lanaceous species, the potato genome sequence database 
were extensively searched using the tomato Phtl genes 
as queries, leading to the identification of a total of 10 
distinct genes as putative potato Phtl genes (Additional 
file 4). Sequence comparison of the potato Phtl genes 
revealed similar amino acid sizes and high sequence 
identities to their corresponding orthologues from to- 
mato (Additional files 4 and 5). It should be noted that 
there also exist two other sequences in the potato gen- 
ome that showed substantial homology to the plant Phtl 
genes, but may be pseudogenes (named as StPTxl and 
StPTx2, respectively), due to the presence of some non- 
sense mutations and the inclusion of some indels (inser- 
tions and deletions) within their putative coding regions. 

Similar to the high sequence identity between the tomato 
and potato Phtl members, high conservation of chromo- 
somal organization of the Phtl homologues from the two 
solanaceous species could also be observed. Figure 2 shows 
the localizations of Phtl genes on the tomato and potato 
chromosomes. It was revealed that the distribution of the 
tomato and potato Phtl genes were obviously uneven, and 
concentrated on only three (3, 6, and 9) chromosomes of 
the two plants. In addition, the supposed three pseudogenes 
{LePTx, StPTxl and StPTx2) were restrictedly assigned on 
the chromosome 9 of the two plants. Interestingly, except 
some individual members, such as LePTS on tomato 
chromosome 6 and StPT9 on potato chromosome 9, most 
of the other Phtl genes/pseudogenes on the corresponding 
chromosomes were distributed in clusters with very short 
physical distance (Additional file 1, Figure 2), suggesting 
that these clustered genes may be produced from in- 
dependent tandem duplications during the evolution of 
Solanaceae Phtl gene family. 

In addition to the Phtl gene themselves, the potential 
genes surrounding each of the Phtl members were also 
carefully surveyed, resulting in the identification of sev- 
eral putative genes exhibiting substantial homology to 
plant Pht3 and Pht4 family genes on the corresponding 
chromosomes (Figure 2). By comparing the locations of 
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Figure 1 Predicted amino acid sequences of the eight tomato Phtl genes LePTl to 8. Sequence alignment analysis was carried out using 
multiple alignment algorithm wrapped within the DNAMAN 7.0 program (http://www.lynnon.com/). Identical amino acids are shaded and gaps 
are indicated by dots. The consensus sites for phosphorylation by protein kinase C and casein kinase II are shown by two arrowheads with red colour 
and the conserved N-glycosylation residue is shown by a green arrowhead. The characteristic Phtl signature was underlined. The transmembrane domains 
(broken underline) were predicted by the Toppred algorithm (httpZ/bioweb.pasteur.fr/seqanal/interfaces/toppred.html). 
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Figure 2 Distribution of Phtl genes [LePTl to 8 and StPTl to 10) on the tomato (T) and potato (P) chromosomes. Chromosome numbers 
are shown at the top of each bar. The arrows next to the gene names indicate the direction of transcription. LePTx, StPTxl and StPTx2 are putative 
Phtl pseudogenes residual in the tomato and potato genomes. The genes from other two families, Pht3 and Pht4, encoding putative Pi transporter or 
carrier with no homology to the Phtl proteins, were also labeled on the corresponding chromosomes of the two plants. 
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these genes, such as LePht3;4 and StPht3;4, we con- 
firmed the existence of two segmental inversions associ- 
ated with the long arms of tomato/potato chromosomes 
6 and 9, which resulted in the inverted linear order of 
the orthologous pairs of PT4/PT5 and PT1/PT3 on the 
corresponding chromosomal regions between the two 
species (Figure 2). 

Phylogenetic analysis of Phtl gene family in tomato and 
other plant species 

In order to perform a comprehensive analysis of evolu- 
tionary relationships among Phtl genes between tomato 
and other plant species, including the eight Phtl proteins 
of tomato, a total of 90 plant Phtl protein sequences, 
representing 11 species from four plant families, Gramineae, 
Brassicaceae, Leguminosae and Solanaceae, were aligned 
and used to construct an unrooted phylogenetic tree. As 
shown in Figure 3, except AtPT6 from Arabidopsis and 
HvPT8, OsPTlS and ZmPTS from each of the three 
graminaceous species, barley, rice and maize, the other plant 



Phtl proteins in the Neighbour-Joining tree were well clus- 
tered into four distinct groups, consisting of one dicot- 
specific group (I), one monocot-specific group (II) and two 
mixed groups (III and IV, respectively) with members from 
both dicots and monocots. 

The Group I harbors the proteins exclusively from the 
dicotyledonous species, in which they are subgrouped by 
phylogeny, and could be further classified into three sub- 
groups (named as subgroup a, b and c, respectively). In 
addition, except the subgroup c, which only includes 
one member from each of the three species, Arabidopsis, 
Medicago and Lotus japonicus, both of the other two 
subgroups, a and b, contain multiple Phtl members 
from the three plant families, Leguminosae, Brassicaceae 
and Solanaceae. For tomato, six of the eight Phtl trans- 
porters (LePTl to 3 and LePT6 to 8) fall into two of the 
three subgroups. Within subgroup a, the two tomato 
members, LePTl and LePT3, group together with their 
orthologous pairs from potato, eggplant and tobacco, to 
the exclusion of other two paralogues, LePT7 and 




Figure 3 Phylogenetic analysis of tomato Phtl genes and other plant Phtl homologs. An unrooted phylogenetic tree of the plant Phtl 
proteins was constructed using the neighbor-joining method with MEGA 5.0 program. Transporters and corresponding plant species are: tomato, 
LePTl to 8 [30,42], this study; potato, StPTl to 10 [30,43,44], this study; tobacco, NtPTl to 5 [45,46]; eggplant, SmPTl to 5 [AS A6\] Arabidopsis 
tholiono, AtPTl to 9 [47]; Medicogo truncotulo, MtPTl to 6 [8,48,49]; Lotus joponicus, LjPTl to 4 [31,50]; Soybean, GmPTl to 14 [13]; Rice, OsPTl to 
13 [28]; Barley, HvPTl to 12 [14,51]; Maize, ZmPTl to 6 [18,29]. 
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LePT8, which group together with other three potato 
homologues and forms a cluster with two soybean ho- 
mologues, GmPT6 and GmPT14. Within the dicot sub- 
group b, the two tomato paralogues, LePT2 and LePT6, 
group closely, and cluster together with three potato ho- 
mologues, StPT2, StPT6 and StPT7. The grouping of 
LePT2 and LePT6 was expected as the two paralogues 
contain identical coding sequences. The rest two paralo- 
gues, LePT4 and LePTS, were found to be assigned only 
into Group III, although the Group IV also contains the 
members from both dicots and monocots. Within Group 
III, the two genes, like their two paralogues LePTl and 
LePTS in Group I, also group together with other solan- 
aceous orthologues by forming an independent Solanaceae 
clade consisting of two subclasses. Both of the subclasses 
contain the orthologous pairs of PT4 or PT5 from tomato, 
potato eggplant and tobacco, suggesting that the duplica- 
tion events associated with the arising of PT4 and PT5, as 
well as PTl and PT3 in tomato and other solanaceous spe- 
cies, occurred before the speciation of Solanaceae lineages. 
Additionally, most of the Phtl proteins in the Group III, in- 
cluding the PT4 and PT5 orthologous pairs, have been ex- 
perimentally evidenced to be strongly induced in the roots 
colonized by arbuscular mycorrhizal (AM) fungi [8]. More- 
over, Phtl members from Arabidopsis, of which the roots 
are unable to form AM symbiosis, are all absent from 
Group III. Interestingly, although the Group IV contains 
much fewer members as compared with the other three 
Groups, there exist two members from each of three spe- 
cies, ArabidopsiSy rice and soybean, but no homologues 
from any of the solanaceous species in the Group IV, sug- 
gesting that the corresponding orthologues from solan- 
aceous lineages have lost after Solanaceae separation with 
Brassicaceae and Leguminosae. 

Expression analysis of the tomato Pht1 genes in different 
tissues under low-P condition 

In this study, for gaining better understanding of the 
possible functions of specific Phtl gene in tomato, the 
tissue-specific expression patterns of each tomato Phtl 
gene were examined in various tissues, including roots, 
stems, young leaves, flowers, as well as fruits at young 
and ripe stages using Real-time RT-PCR. The quantita- 
tive data showed that except LePT4 and LePTS, of which 
the transcripts were not detectable in all tissues exam- 
ined, the transcripts of other Phtl paralogues were all 
detectable in a certain tissues and showed distinct but 
partially overlapping expression profiles (Figure 4). 

LePTl was expressed in all tissues examined, and its 
transcripts were detected abundantly in roots and leaves, 
and to a lesser extent in stems and flowers, as well as in 
fruits. The transcripts of LePTl in green fruits were four 
times more than those in ripe fruits. In contrast to the 
ubiquitous expression profiles of LePTl, the expression 



of LePT2 showed relatively distinct tissue-specific pro- 
files, with its transcripts intensively in roots and ex- 
tremely faintly in some of other tissues, such as in green 
and ripe fruits. The expression patterns of the LePTS 
and LePTS were a little similar, as both of the two genes 
were expressed very weakly in all tissues. Even so, the 
highest transcript level for LePTS was detected in ripe 
fruits, and was about ten times more than that in green 
fruits. LePT6, the closest fellow of LePT2 in phylogeny, 
was also dominantly expressed in roots, but with only 
one-third of the expression level of LePT2 in the root 
tissues. Additionally, very weak transcript levels of this 
gene were also detectable in stems and leaves. LePT7 
was also ubiquitously expressed in all tissues, and had a 
very similar expression tendency, but significant lower 
expression levels in all tissues as compared to its paralo- 
gue, LePTl (Figure 4). The differential but overlapping 
expression of the Phtl genes well mirrors the evolution- 
ary conservation and functional divergence of Phtl 
transporters in tomato plants. 

Expression analysis of tomato Phtl genes in response to 
AMF colonization under low and high Pi supply 
conditions 

Since the expression of some Phtl genes in tomato and 
also in other plant species, have been characterized to be 
AM-inducible and Pi-responsive [21,52], the relative ex- 
pression levels of each tomato Phtl member were thus 
further determined in roots and leaves in response to 
AM Fungi (Glomus intraradices) colonization under low 
(50 (iM) and high (1 mM) Pi supply condition. As shown 
in Figure 5, colonization of AM fungi increased not only 
the biomass, but also the P concentration of the tomato 
plants under the low Pi supply condition; however, no 
significant difference of both the biomass and P concen- 
tration could be observed between the colonized and the 
noncolonized plants under the high Pi supply condition. 

qRT-PCR analysis revealed that except the three para- 
logues, LePTS, LePT4 and LePTS, of which the transcripts 
were strongly enhanced or specifically activated only in 
the inoculated roots under the low Pi supply condition, 
and LePTS, of which the transcripts were not detectable in 
both tissues under any treatments, the expression of the 
other four paralogues, LePTl, LePT2, LePT6 and LePT7, 
were significantly repressed under the high Pi supply con- 
dition (Figure 6). Such down-regulated cases occurred 
more conspicuously upon the two paralogues, LePT2 and 
LePT7, as their transcripts in both the root and leaf tissues 
were drastically decreased {LePT2) or even completely ab- 
sent (LePTT) under high Pi conditions regardless of with 
or without AM colonization. In addition, very significant 
decrease of the transcript abundance of the four paralo- 
gues was also detected in both the roots and leaves of the 
colonized tomato plants as compared to those non- 
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Figure 4 Tissue-specific expression analysis of tomato Phtl genes. The RNA were prepared from different tissues, including roots (R), stems 
(S), young leaves (L), flowers (FL), as well as fruits at green (GF) and ripe (RF) stages. The relative expression levels of each of the tomato Phtl 
genes were indicated as percentage of the constitutive Actin expression activity. Each bar was the mean of three biological replications with 
standard error. 



colonized controls under low Pi supply condition (Figure 6). 
The remarkable down regulation of these four members in 
response to high-P supply and AMF-colonization might be 
partially caused by the significant increase of P concentra- 
tion in such treated tomato plants (Figure 5). Interestingly, 
although LePT2 and LePT6 were considered to be the clos- 
est related genes in tomato Phtl family due to their identi- 
cal coding sequence, the down regulation of LePT6 in roots 
in response to AM symbiosis under low Pi condition was 
much moderate than that of LePT2, Such discrepancy in 
expression levels strongly suggests that the regulatory com- 
ponents controlling the activation or suppression of LePT2 
and LePT6 have divergent after the two paralogues pro- 
duced from a relatively recent duplication event. 

The specialized expression profiles of the tomato Phtl 
genes in response to AM symbiosis or different Pi status 
prompted us to investigate their promoter regions. As 
shown in Figure 7A, the number and localization of the 
two Pi-regulated (PIBS and W-box) [53,54] and one 
AM-responsive elements (MYCS) [45,55] differ widely in 
the promoter regions of these eight Phtl genes, even 
though the coding sequence and expression profiles of 



some paralogues are highly conserved. However, similar 
as the AM-induced Phtl genes in other dicot species, 
the MYCS motif was found to be present exclusively in 
the putative promoter regions of the three AM-activated 
Phtl paralogues, LePT3, LePT4 and LePTS, and were lo- 
cated very closely to the Pi-regulated PIBS element [45]. 
Histochemical staining analysis further revealed that the 
LePTS and LePTS promoter regions {pLePT3_225o and 
pLePTS^4yi) containing the two elements, MYCS and 
PIBS, were sufficient to direct p-glucuronidase (GUS) 
expression specifically in the mycorrhizal roots and were 
limited to distinct cells harboring AM fungal structures 
(arbuscules or intracellular hyphae) (Figure 7B), similar 
to the cellular distributions of their paralogue LePT4 
and other AM-inducible Phtl homologues from various 
other plant species reported previously [30,55-57]. 

Discussion 

In recent studies, benefiting from the availability of 
whole genome sequence of model plants, dozens of 
genes belonging to the Phtl family that encode putative 
high-affinity Pi transporters have been identified from 
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various plant species using the comparative genome ap- 
proaches. Tomato, a model species from the Solanaceae 
family, has been historically characterized to have at 
least five Phtl genes [58]. Our present study, through 
extensive searches of available databases, led to the iden- 
tification of a total of eight putative Phtl genes in the to- 
mato genome. As this is the first genome-wide analysis 
of the Phtl gene family in any solanaceous species, the 
investigation of chromosomal organization, evolutionary 
relationships, as well as expression patterns of the to- 
mato Phtl genes in this study is of great significance and 
would offer a basis for better understanding the evolu- 
tionary mechanisms underlying the expansion, conserva- 
tion and functional divergence of the Phtl genes in the 
whole Solanaceae family. 

Evolutionary expansion of the tomato Phtl genes 

Multigene families, in a general way, could arise through 
tandem duplications, resulting in a clustered occurrence, 
or through genome/segmental duplications, resulting in 
a discrete distribution of family members [59]. As most 
of the tomato Phtl genes were assigned in clusters (such 
as PT1/PT3, PT2/PT6 and PT4/PT5), with not only very 



close physical localization (Figure 2), but also very high 
levels of sequence identity (Table 1, Additional file 5), it 
is strongly suggestive of that tandem duplications might 
be the major contributors to the expansion of the to- 
mato Phtl family. Additionally, since most of the tomato 
Phtl members group together with their orthologues 
from other solanaceous species, such as potato, eggplant 
and tobacco by forming independent solanaceous clades 
to the exclusion of other dicot homologues (Figure 3), 
indicates that the duplications associated with the arising 
of the coupled paralogues such as PT1/PT3 and PT4/ 
PT5 in solanaceous species, occurred before the speci- 
ation of solanaceous lineages from a common ancestor. 
Intriguingly, in viewing of the localization of tomato 
LePT2/6 and potato StPT2/6/7 on their corresponding 
chromosomes, it is tempting to make a tendentious con- 
clusion that the duplication giving rise of the tomato 
LePT2 and LePT6 probably occurred before the split of 
tomato and potato. However, the distribution of the 
LePT2 and LePT6 in the terminal subclades of the 
phylogenetic tree, and the identical coding sequence 
shared by them well reflected that the two paralogues 
were produced from the more recent duplication events 
that occurred within the tomato lineage postdating it 
split from a common ancestor shared by potato. 

It has been recently documented that the Solanum 
lineage genome has undergone two rounds of consecutive 
whole-genome triplication events, one that was ancient and 
shared with most dicot plant families, and one that was 
more recent and occurred before the divergence of tomato 
and potato lineages [34], which led to the hypothesis that 
segmental duplications produced by genome polyploidy 
may also exert important impact on the expansion of the 
Phtl family. The loci of LePT2/6 on chromosome 3 and 
LePT4/5 on chromosome 6 flanked respectively by two 
paralogues, LePht3;3 and LePht3;4 (with non-homology to 
the Phtl genes) (Figure 2) strongly suggests that the arising 
of the two pairs, PT2/PT6 and PT4/PT5 might originate 
from a segmental duplication, followed by two independent 
tandom duplications, which eventually resulted in the fix- 
ation of the two couples of Phtl members on the chromo- 
somes 3 and 6. As the PT4 and PT5 paralogues cluster 
together with other members from both dicots and mono- 
cots in the phylogenetic tree (Figure 3), indicating that the 
segmental duplication yielding the precursors of the two 
couples, PT2/PT6 and PT4/PT5, occurred before the diver- 
gence of monocots and dicots. With regard to the other 
four paralogues, the clustered LePTl and LePT3, and the 
individual LePT7 and LePT8, they may be the result of sev- 
eral relatively recent segmental or single-gene duplication 
events that occurred before the speciation of tomato and 
potato lineages, and followed by at least one independent 
tandom duplication event (producing the two paralogues, 
LePTl and LePT3), It has been well documented that 
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genome polyploidizations are commonly accompanied by 
massive chromosomal rearrangements [60]. In our study, 
by comparing the phylogenetic tree and the chromosomal 
distribution of the tomato and potato Phtl genes, two seg- 
mental inversions leading to the inconsistent linear orders 
of PT1/PT3 and PT4/PT5 were identified between the two 
solanaceous genomes, well supporting the very recent find- 
ings that at least nine large and several smaller inversions 
exist between the tomato and potato lineages [41]. 

Functional conservation and divergence of the tomato 
Phtl gene family 

It has been generally accepted that gene duplication followed 
by functional differentiation has performed a pivotal role in 
driving evolutionary novelty that allow plants to increase fit- 
ness to new environments [47]. To data, as the lack of 
genome-wide survey of Phtl genes in any solanaceous spe- 
cies, there is no systematic analysis of tissue-specific expres- 
sion patterns for the tomato Phtl family so far. In our 



present work, we revealed that differential but partial over- 
lapping expression of the Phtl genes occurred in tomato, as 
did the members of this family in several other plant species, 
such as ArabidopsiSy rice and soybean [13,28,43]. The spe- 
cialized expression of these genes well mirrors the evolution- 
ary divergence of regulatory elements that are required for 
controlling Pi uptal<e and mobilization within/across particu- 
lar tissues or cells during tomato plant growth. 

Earlier results, based on the study of tissue-specific ex- 
pression and cellular distribution of Phtl genes in several 
different plant families revealed that many of the Phtl 
genes are expressed dominantly in roots, especially in root 
epidermis and root hairs, in response to P deprivation, 
suggesting a potential role of these genes in Pi capture and 
uptake [12,19]. The transcription data obtained in this 
study indeed provide direct evidence for strong expression 
of most of the tomato Phtl genes in the roots under low 
Pi supply condition (Figure 6). It was shown that al- 
though the transcripts of LePTl in the roots and leaves 
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significantly decreased in response to Pi sufficiency, a low 
level of constitutive expression could be detected through- 
out the plant, consistent with the expression patterns of 
its orthologue, StPTl, in potato [46], suggesting that 
LePTl and its orthologues may be involved in not only up- 
take of Pi from soil solution but also redistribution of Pi 
within plants. LePT2 has been previously documented to 
be expressed exclusively in P-depleted roots. However, in 
the present work, a relatively weak but still observable 
transcription level could be detected in the roots irrigated 
with high Pi (1 mM) solution, similar results could also be 
observed from our previous studies on the LePT2 ortholo- 
gues in other three solanaceous species, eggplant, pepper 
and tobacco [61]. In addition, very slight levels of the LePT2 



transcripts were also detectable in stems, flowers and fruits 
at green and ripe stages under low Pi supply condition. Al- 
though LePTl shares its coding sequence identical to its 
paralogue, LePTS, the transcript abundance of the two 
members were observably different whether under low Pi 
supply condition or in response to AM symbiosis. Such dis- 
crepancy between the two close paralogues may be caused 
by the inconsistent distributions of Pi-responsive elements, 
such as PIBS and W-box, in their promoters (Figure 7A) 
[53,54,62-64]. Even so, the identical protein activity and high 
degree of overlapping expression strongly implies the pres- 
ence of functional redundancy between the two members. 

With regard to the three AM-activated Phtl paralo- 
gues, LePTS, LePT4 and LePTS, as their transcripts 
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could be induced abundantly only in the inoculated 
roots under low Pi, but not high Pi condition, similar as 
the expression of their orthologues from tobacco, pep- 
per, and eggplant, suggesting that the AM-activated ex- 
pression of LePTS'S, like their solanaceous orthologues, 
might also be regulated by at least two conserved c/5-ele- 
ments, PIBS and MYCS [45]. Comparative screening of 
the putative promoter regions of the tomato Phtl genes 
indeed led to the identification of the two motifs that 
were closely and exclusively localized in the promoters 
of the three AM- activated Phtl genes (Figure 7A). The 
relatively high proportion of AM-induced Phtl genes 
present/retain in solanaceous species compared to that 
in many of other mycorrhizal plants, such as rice, to a cer- 
tain extent, reflects the importance of these genes in regula- 
tion of symbiosis during the Solanaceae evolution. In 
addition, the strong expression levels of these symbiosis- 
activated Phtl genes in contrast to the remarkable down- 
regulation of the Pi transporters responsible for direct Pi 
uptake in mycorrhizal roots also provided strong evidence 
to support the earlier findings that symbiotic uptake path- 
way would contribute the majority of the accumulated Pi 
received by mycorrhizal tomato and other plants [65-67]. It 
is worth noting that although LePT4 was observed to be 
the gene that had the highest expression level than the 
other paralogues in tomato mycorrhizal roots under low Pi 
condition (Figure 6), mutation of the LePT4 expression in 
tomato virtually unaffected the establishment of AM sym- 
biosis [30,40], which seems contradictory to the recent 
findings that knock down/out of the AM-specific or up- 
regulated Phtl genes in Lotus japonicus, Medicago and 
Astragalus sinicus significantly impaired both the deve- 
lopment of AM interaction and symbiotic Pi uptake 
[31,55,68]. The absence of AM-associated phenotypes in 
the tomato lept4 knock out mutant was thus suggested to 
be the genetic redundancy within the tomato Phtl gene 
family [30,69]. This explanation might be reasonable as a 
result of that LePT4 and its AM-induced paralogue, LePTS, 
were considered to be diverged from a common precursor 
through tandom duplications (Figure 2), and thus it is un- 
surprising that LePTS would share a similar physiological 
role with LePT4 in tomato. Interestingly, although there 
exist two AM-specific Phtl genes, the strongly activated 
OsPTll and poorly induced OsPTlS, in the rice Phtl fam- 
ily, no fiinctional redundancy was observed between the 
two paralogues, as silencing whichever of the two paralo- 
gues caused a significant repression of AM symbiosis [70]. 
Since the rice OsPTll and OsPT13 are distributed rela- 
tively distantly in phylogenetic tree, similar as the phylogen- 
etic relationships between solanaceous PT3 and PT4/PT5 
paralogues, it is tempting to raise a speculation that no 
fijinctional redundancy might be existent between the two 
groups (PT3 and PT4/PT5) of Phtl transporters in regula- 
tion of the development of AM symbiosis. 



Interestingly, within some other plant species, such as in 
soybean, there exists specialized Phtl member(s) that is 
(are) mainly expressed in sink tissues, such as in flowers, im- 
plying possible fijinctions of these genes in Pi import from 
source to sink [13,71]. However, in our present study, we 
did not observe any member in the tomato Phtl family that 
is dominantly expressed in the flower or fruit tissue. The 
relatively high transcription levels of LePTl and LePT7 in 
these sink tissues compared to the other six paralogues led 
to the suggestion that the two members, especially the 
LePTl, might have been evolved to meet the requirement of 
transporting Pi from the source to the sink organs or cells. 
One of the potential evolutionary fates of gene duplications 
is considered to silence (nonfunctionalization) one of the 
duplicate copies [72]. The inactivation of LePTS in all the 
tissues led to the suggestion that LePTS might be on the 
way to become a pseudogene. The presence of several frag- 
ments that showed substantial homology to plant Phtl 
genes but with incomplete coding regions in the tomato and 
potato genomes also well supports the theory that genome 
polyploidizations as well as the following gene loss (diploidi- 
zation) are common characters of plant genomes [73,74]. 

Conclusions 

Taken together, this study provided the first comprehen- 
sive analysis of the chromosomal organization, phylo- 
genetic evolution and tissue-specific expression patterns 
for each member of the Phtl family in tomato. The re- 
sults presented here could offer a useful basis for future 
research work on better understanding the mechanisms 
underlying the evolutionary regulation of Phtl genes in 
response to Pi deficiency and AM symbiosis during to- 
mato growth. The high conservation not only in the 
coding sequence, but also in the chromosomal distribu- 
tion between the tomato and potato Phtl orthologues 
could also lend strong evidence to support the further 
comparative genomics analysis across the whole Solanaceae 
family. However, we also realize that although Phtl family 
have been commonly considered to be a high-affinity Pi 
transporter family, an increasing number of members of 
this family in other plant species have been demonstrated 
exhibiting a low or even dual affinity for Pi uptake in heter- 
ologous yeast or oocyte expression system [75,76]. There- 
fore, more information, especially the transport kinetics, 
cellular distributions and physiological phenotypes of knock- 
out/knockdown mutants is needed in the near fixture to de- 
termine more precise fijinctional roles for each of the tomato 
Phtl genes. 

Methods 

Plant material and growth conditions 

Tomato {Solarium lycopersicum cv. Micro-Tom), was 
used in this study. The seeds were surface-stenlized, ger- 
minated and maintained in tissue culture for three weeks 
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using 1/2 MS medium supplemented with 1.5% sucrose. 
The aseptic plantlets were then transferred to pot cul- 
ture for either with high-/low-P treatment or inoculation 
with AM fungi. 

In pot culture, two plantlets were transplanted to a 3 dm^ 
plastic pot filled with sterilized sand. A sand-based inocu- 
lum containing Glomus intraradices {Gi) was used for in- 
oculation. Each plants was inoculated with 5 g inoculum or 
autoclaved inoculum around the roots. The plants were 
grown in a growth room with a 14-h light period (28-30°C) 
and a 10-h dark period (18-20°C). The irrigating solution 
contained the following nutrients: 1 mM NH4NO3, 2 mM 
KNO3, 0.5 mM Ca(N03)2, 0.25 mM CaCb, 0.5 mM 
MgS04, 20 (iM Fe-EDTA, 9 \iM MnCb, 46 (iM H3BO3, 
8 (iM ZnS04, 3 (iM CUSO4, 0.03 (iM (NH4)2Mo04, and 
1 mM (high-P treatment) or 0.05 mM Pi (low-P treatment) 
NaH2P04. The experiment comprised 4 replicates for each 
treatment. After inoculation for six weeks, the plants were 
either harvested for collecting the root, stem and young leaf 
samples or for continuing growing for the later collection 
of flower and fruit samples at young and ripe stages. The 
collected samples were immediately frozen in liquid nitro- 
gen and stored at -80°C for subsequent RNA isolation. 

Identification of Phtl genes in tomato and potato genomes 

Members of Phtl gene family in the tomato genome 
were identified using the BLASTN and TBLASTN algo- 
rithm wrapped in the BLAST 2.2.27+ applications. To 
identify the potential Phtl genes in tomato genome, cod- 
ing sequence (cds) and deduced protein sequences of 
the Arabidopsis and rice Phtl genes were queried re- 
spectively in the tomato genomic sequence database 
downloaded from the Solanaceae Genomics Network 
(www.sgn.cornell.edu). Sequences with a query over 50% 
and e-value less than -10 were taken as the Phtl candi- 
dates. AU the obtained sequences were submitted to 
NCBI (http://www.ncbi.nlm.nih.gov/) and Pfam database 
(www.sanger.ac.uk/Software/Pfam/search.shtml) for fur- 
ther confirmative analysis. For chromosomal localization 
analysis, the tomato Phtl candidates were further used 
as queries for BLASTN searches against the SGN To- 
mato Whole Genome Scaffolds data (2.40) (http://solge- 
nomics.net/organism/Solanum_lycopersicum/genome). 

For identification of the potential homologues of tomato 
Phtl genes fi:om potato genome, the potato genome sequence 
data downloaded fi:om the SGN database (http://solgenomics. 
net/organism/Solanum_tuberosum/genome) were also exten- 
sively searched using the tomato Phtl genes as queries. The 
naming of the potato Phtl genes were partially based on their 
phylogenetic relationships with the tomato homologues. 

Phylogenetic analysis 

The sequence data used in this study were coUected 
using a query search in the NCBI database using the 



known Phtl family gene sequences from Arabidopsis 
and rice. Multiple sequence alignments were performed 
using the program ClustalX (version 1.8) with default 
gap penalties. An un-rooted phylogenetic tree was gen- 
erated using the deduced amino acid sequences of Phtl 
genes by neighbor-joining algorithms wrapped in the 
MEGA 5.1 phylogeny program (www.megasoftware.net). 
Bootstrap analysis was carried out with 1,000 replicates. 

RNA extraction and real-time RT-PCR analysis 

Total RNA was isolated from 100 mg of various tissue 
samples, including roots, stems, young leaves and 
flowers, using the guanidine thiocyanate extraction 
method with Trizol reagent (Invitrogen) and from fruit 
samples using CTAB-sour phenol extraction method as 
described by Chang [77]. After extraction, the RNA sam- 
ples were treated with DNase I (TaKaRa) to eliminate 
the trace contaminants of genomic DNA. For conduct- 
ing reverse transcription (RT) PCR analysis, approxi- 
mately two micrograms of total RNA from each sample 
was used to synthesize first-strand cDNA using a reverse 
transcription kit (TaKaRa), and the synthesized cDNAs 
were used as templates in the following PCR reactions. 

Real-time RT-PCR analysis was performed to relatively 
quantify the expression levels of tomato Phtl genes in 
different tissues or in roots and leaves in response to 
mycorrhizal colonization at high and low P status. The 
reaction was conducted on the Applied Biosystems 
(ABI) Plus Real-Time PCR System using the SYBER pre- 
mix ExTaq kit (TaKaRa). Relative quantification of the 
transcripts for each tomato Phtl gene was standardized 
to the expression level of the tomato constitutive Actin 
gene, calculated by the formula Y = 10'^"^^^^^^ x 100% 
(ACt is the differences of cycle threshold value between 
the target Phtl gene and the control Actin products) 
[47,78]. The specificity of primer sets designed for the 
qRT-PCR (Additional file 6) was confirmed by sequen- 
cing after the PCR reaction. 

Histochemical GUS staining and detection of mycorrhizal 
fungal colonization 

A 1250-bp LePTS promoter fragment and a 471 -bp LePTS 
promoter fragment immediately upstream of the transla- 
tion initiator ATG were amplified and cloned into binary 
vector pBI12, respectively, to replace the CaMV35S pro- 
moter in front of the |3 -glucuronidase (GUS) gene. The 
resulting two constructs were designated as pLePT3_225o 
and pLePTS_4yi, respectively, and introduced into 
Agrobacterium tumefaciens strain EHA105 for gen- 
etic transformation. 

Histochemical GUS staining of the fresh transgenic roots 
was performed as described previously [44]. For visuali- 
zation of fiangal structures, the Magenta-GUS stained root 
segments were treated with 10% KOH solution heated to 
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90°C for 1 h, and then neutralized with 1% HCl (v/v) solu- 
tion for 5 min. The root materials were then counter- 
stained with 0.3% trypan blue solution for 2 h at 90°C. The 
co-localization of Magenta-GUS and the trypan blue stain- 
ing were indicated by purple color. The stained materials 
were rinsed in 50% glycerol and photographed by a stereo- 
microscope with a color CCD camera (Olympus). 

Additional files 
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